Isotherms clustering in cosmic microwave background 
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Isotherms clustering in cosmic microwave background (CMB) has been studied using the 3-year 
WMAP data on cosmic microwave background radiation. It is shown that the isotherms clustering 
could be produced by the baryon-photon fluid turbulence in the last scattering surface. The Taylor- 
microscale Reynolds number of the turbulence is estimated directly from the CMB data as Rex ~ 

PACS numbers: 98.70.Vc, 98.80.-k, 98.65.Dx 
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INTRODUCTION 

Luminous matter in the observable universe is strongly 
clustered (see, for a review jlj). Main reasons for this 
apparent clustering are gravitational forces and an initial 
non-uniformity of the baryon-photon matter just before 
the recombination time. An information about this 
non-uniformity we can infer from the Cosmic Microwave 
Background (CMB) radiation maps. Although these 
maps indicate nearly Gaussian distribution of the CMB 
fluctuations, certain clustering may be observed already 
in these maps even after they have monopole and dipole 
set to zero. The clustering of the isotherms in the 
last scattering surface then can be used as an initial 
condition for the gravitational simulations. One of 
the purposes of this paper is to describe the clustering 
quantitatively. Another related purpose is to find out 
underlying primordial physical processes, which may 
cause the isotherms clustering in the CMB maps. It is 
known that rotational velocity perturbations in 

the primordial baryon-photon fluid can produce angular 
scale anisotropics in CMB radiation through the Doppler 
effect. The conclusions are relevant to arcminute scales 
(for large scales see 4]). On the other hand, it is recently 
shown in 5j that clustering is an intrinsic property of 
the turbulent (rotational) velocity field. Therefore, one 
should check whether the CMB isotherms clustering can 
be related to turbulent motion of the baryon-photon 
fluid just before recombination time (an observational 
indication of turbulent motion in the baryon-photon 
fluid just before recombination time is given in [g, see 
also Moreover, we will estimate the Taylor- 

microscale Reynolds number {Rex, see Appendix A) of 
this motion, that is Re\ ~ 10^ (let us recall that critical 
value is Rex ~ 40 ). As far as we know this is first 
estimate of the Reynolds number obtained directly from 
a CMB- map (we have used a cleaned 3-year WMAP 0|). 



fluctuations are dominated by acoustic perturbations of 
the primordial velocity field. The velocity field can be 
represented by following way 



where the acoustic component is a potential one Ua = 
and Ur is rotational component of the primordial 
velocity field. To exclude the acoustic component one 
can take operation curl u = curl Ur, that is vorticity. 
Average energy dissipation in turbulent velocity field is 



determined just by its rotational component 
{e) v{{curl uf) 
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where (...) means statistical average and v is kinematic 
viscosity. Since it is expected that high frequency events 
in the velocity field should provide the most significant 
contribution to the turbulent dissipation and, especially, 
to its high order moments one can also expect that clus- 
tering in turbulent velocity field should be intimately re- 
lated to fluctuations of the ener gy d issipation (so-called 
intermittency phenomenon 

Il31l3)- Taking into ac- 
count that turbulent velocity is nearly Gaussian 

[13 (see 

also Appendix B) the clustering phenomenon in the tur- 
bulent velocity field should be also a Gaussian one. On 
the other hand, the intermittency phenomenon in fluid 
turbulence, usually associated with fine non-Gaussian 
properties of the velocity field Thus, simultane- 

ous consideration of the cluster and intermittency char- 
acteristics can provide an additional valuable information 
about physical origin of the isotherms clustering in the 
CMB maps and about the fine non-Gaussianity of the 
maps. Therefore, the paper starts with an introduction 
to this subject. 



CLUSTERING AND INTERMITTENCY 



To find physical (dynamical) origin of the isotherms 
clustering in the CMB we will compare certain statis- 
tical properties of the CMB maps (and their Gaussian 
simulations) with corresponding properties of the fiuid 
turbulence. It is well known that energetically the CMB 



For turbulent flows the so-called Taylor hypothe- 
sis is generally used to interpret the turbulent data 
18] . This hypothesis states that the intrinsic 
time dependence of the wavefields can be ignored when 
the turbulence is convected past the probes at nearly 
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FIG. 1: Standard deviation of the running density fluctua- 
tions against r for for the velocity fluctuations measured in 
the wind-tunnel (Fig. la) and for their Gaussian spectral 
simulation (Fig. lb). 
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FIG. 2: As in Fig. 1 but for atmospheric surface layer at 
Re\ = 19000 (Fig. 2a) and for its Gaussian spectral simula- 
tion (Fig. 2b). 



constant speed. With this hypothesis, the temporal dy- 
namics should reflect the spatial one, i.e. the fluctuating 
velocity field measured by a given probe as a function of 
time, u{t) is the same as the velocity u{x/ {u)) where {u) 
is the mean velocity and x is the distance to a position 
"upstream" where the velocity is measured at t = 0. 

Let us count the number of 'zero'-crossing points of 
the signal in a time interval r and consider their running 



density n^-. Let us denote fluctuations of the running 
density as Jrir = n-r — {tIt), where the brackets mean 
the average over long times. We are interested in scaling 
variation of the standard deviation of the running density 
fluctuations {Sn'^)^^''^ with r 
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For white noise it can be derived analytically |l9j,|20j 
that a„ = 1/2 (see also 15|). 

In Figure la and 2a we show calculations of the stan- 
dard deviation for a turbulent velocity signal obtained 
in the wind-tunnel experiment |2ll | and for a velocity 
signal obtained in an atmospheric experiment |22l| ) re- 
spectively. The straight lines are drawn in the figures 
to indicate scaling (1). One can see two scaling inter- 
vals in Fig. la. The left scaling interval covers both 
dissipative and inertial ranges, while the right scaling in- 
terval covers scales larger then the integral scale of the 
flow. While the right scaling interval is rather trivial 
(with Un = 1/2, i.e. without clustering), the scaling in 
the left interval (with a„ < 1/2) indicates clustering of 
the high frequency fluctuations. The cluster-exponent 
Un decreases with increase of Re\, that means increas- 
ing of the clustering (as it was expected from qualitative 
observations). 

To describe intermittency the standard deviation of 
fluctuations of the energy dissipation is used [i^l (see 
also a discussion in l23 |) 



2\l/2 



where 



or in the terms of the discrete time series ^3 



(2) 



(3a) 



(36) 



where Auj = Wj+i — Uj. 

Figure 3 shows an example of the scaling (2) providing 
the intermittency exponent ai for the wind-tunnel data 
(see |2^ for more details). Events with high concentra- 
tion of the zero-crossing points in the intermittent tur- 
bulent velocity signal provide main contribution to the 
(fc^) due to high concentration of the statistically sig- 
nificant local maximums in these events |25j . Therefore, 
one can make use of statistical version of the theorem 
'about mean' in order to estimate the standard deviation 
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2}2)l/2(^„2^1/2 



(4) 



where {max : u^}t is maximum of w(i)^ in interval of 
length T. Value of {S{max : u^}^)^/^ should grow with 
Re\ (due to increasing of (u^)) and, statistically, with 
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FIG. 3: Intermittency exponent as slope of the straight hne FIG. 5: The cluster-exponent an (circles) and intermittency 
for the longitudinal velocity signal obtained in the wind tunnel exponent Ui (squares) against 1/ In Rex for velocity signals at 
experiment i2ll at Rex ~ 1029. different values of Reynolds number (200 < Rex < 20000). 
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FIG. 4: Statistical ('tail') exponent Om — 0.1 as slope of 
the straight line for the velocity signal obtained in the wind 
tunnel experiment [2H at Rex = 1029. 



the length of interval r (the later growth is expected to 
be self-similar), i.e. 



(J{max:M2}2)i/2 ^C'(i?eA) 



(5) 



where constant C{Re\) is a monotonically increasing 
function of Re\, and the statistical ('tail') exponent 
dm > is independent on Re\. Fig. 4 shows an example 
of scaling (5) providing the statistical ('tail') exponent 
am — 0.10 ± 0.01 for the wind-tunnel data. Substituting 
(5) into (4) and taking into account (1) and (2) we can 
infer a relation between the scaling exponents 
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(6) 



Unlike the exponent a™ the exponent a„ depends on 
Re\. Therefore, we can learn from (6) that just the clus- 
tering of high frequency velocity fluctuations (a^ in (6)) 
is responsible for dependence of the intermittency expo- 
nent ai on the Re\. 

It is naturally consider the exponent functions through 
In R\ : a = a (In Rex) li.lii.lia (see also Appendix A). 
Following to the general idea of Ref. I^^l (see also [^j ) let 
us expend this function for large Rex in a power series: 



Q!(ln Rex) — a{oo) 



ai 



+ 
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lni?eA (InRex)^ 



(7) 



In Figure 5 we show the calculated values of Q!„ (circles) 
and Ui (squares) for velocity signals against 1/ In Re\ for 
200 < Re\ < 20000 (for ai we also added four data points 
taken from the paper , where the data were obtained 
in a mixing layer and in an atmospheric surface layer). 
The straight lines (the best fit) indicate approximation 
with the two first terms of the power series expansion 
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This means, in particular, that 
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In Re\ 







(8) 



(9) 



The closeness of the constants oi ~ 3/2 (the slopes of the 
straight lines in Fig. 5) in the approximations (8) of the 
cluster-exponent q;„ and of the intermittency exponent 
ai confirms the relationship (6) with am — 0.1 (cf. Figs. 
4 and 5). 

Thus one can see from Eq. (6) that entire dependence 
of the intermittency exponent ai on Rex is determined by 
dependence of the cluster-exponent an on Re\. On the 



4 



other hand, the cluster-exponent itself is uniquely deter- 
mined by energy spectrum of the velocity signal at sug- 
gestion that the velocity field is Gaussian. However, dif- 
ferent spectra can produce the same value of the cluster- 
exponent (in particular, the cluster-exponent is invariant 
to variation of the potential component of the velocity 
field, that is significant for next section). 

We use spectral simulation to illustrate that en- 
ergy spectrum uniquely determines cluster-exponent for 
Gaussian-like turbulent signals. In these simulations we 
generate Gaussian stochastic signal with energy spec- 
trum given as data (the spectral data are taken from 
the original velocity signals). Figures lb and 2b show, as 
examples, results of such simulation for the wind-tunnel 
velocity data (Fig. lb) and for the atmospheric surface 
layer (Fig. 2b). Namely, Figures la and 2a show results 
of calculations for the original velocity signals and Figs, 
lb, and 2b show corresponding results obtained for the 
spectral simulations of the signals. This means, in par- 
ticular, that in respect to the cluster-exponent the tur- 
bulent velocity signals exhibit Gaussian properties. Of 
course, it is not the case for the intermittency exponent 
ai. Namely, at ~ 0.5 for the Gaussian spectral simula- 
tions, independently on the spectral shape. 



CMB CLUSTERING AND INTERMITTENCY 

The motion of the scatterers in the last scattering 
surface imprints a temperature fluctuation, AT, on the 
CMB through the Doppler effect 0,(3 



AT(n) 
T 



J g{L)n ■ vb(x) dL 



where n is the direction (the unit vector) on the sky, Vb 
is the velocity field of the baryons evaluated along the 
line of sight, x = Ln, and g is the so-called visibility. It 
should be noted that, in a potential flow, fluctuations per- 
pendicular to the line of sight lack a velocity component 
parallel to the line of sight; consequently, generally there 
is no Doppler effect for the potential flows. The same 
is not true for rotational flows, since the waves that run 
perpendicular to the line of sight have velocities parallel 
to the line of site 

Clustering in the turbulent velocity field in the last 
scattering surface, therefore, can produce clustering of 
the observed CMB isotherms. Thus (see previous sec- 
tion), the areas of strong clustering of the CMB isotherms 
can indicate the areas of the strong vorticity activity in 
last scattering surface. Moreover, using results of pre- 
vious section we can, calculating corresponding cluster- 
and intermittency exponents for the CMB maps, check 
whether the turbulent relationship (6) is also valid for 
the CMB data, with the same value of am ~ 0.1 as for 
the fluid turbulence. In the case of positive answer we 



can then estimate the Taylor-microscale Reynolds num- 
ber Re\ for the primordial turbulence in the last scatter- 
ing surface using Fig. 5. 

We used the 3-year WMAP data cleaned from fore- 
ground contamination by Tegmark's group (the orig- 
inal temperature is given in fiK units). Let us introduce 
a gradient measure for the cosmic microwave radiation 
temperature T fluctuations that is an analogy of the tur- 
bulent energy dissipation measure 0|,[l3 



Xr 



(10), 



where Vr is a subvolume with space-scale r (this is a 
space analogy of Eqs. (3a, b), see also Introduction: 
(e) ~ {{curl v)2)). 
Then, scaling law 



{Sxl 



(11) 



corresponds to the scaling law (2) and the exponent ai is 
the CMB intermittency exponent. Technically, using the 
cosmic microwave pixel data map, we will calculate the 
cosmic microwave radiation temperature gradient mea- 
sure using summation over pixel sets (by random walks 
method) instead of integration over subvolumcs Vr- The 
scaling of type (11) (if exists) will be then written as 



{Sx 



2\l/2 



(12) 



where the metric scale r is replaced by number of the pix- 
els, s, characterizing the size of the summation set (the 
random walk trajectory length). The Xs is a surrogate of 
the real 3D gradient measure Xr- It is believed that such 
surrogates can reproduce quantitative scaling properties 
of the 3D prototypes 0| . Analogously to the previous 
section the cluster analysis can be performed 
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(13) 



Figure 6a shows the c^Mster-scaling (13) for the 
CMB calculated for the 3-year WMAP map cleaned by 
Tegmark's group The straight line (the best fit) is 
drawn to indicate the scaling in log-log scales. Figure 
6b shows the results of analogous calculations produced 
for Gaussian spectral simulation of the 3-year WMAP 
map. As it was expected the Gaussian spectral simula- 
tion produces map with the same cluster properties as 
the original map (cf. previous section). 

The intermittency exponent ai, on the contrary, takes 
significantly different values for the original CMB map 
and for its Gaussian spectral simulation (Figs. (7a, b)). 
One can see, that for the Gaussian map-simulation ai ~ 
0.5 ± 0.02 (as it is expected, see previous section), while 
for the original CMB map ai ~ 0.36 ± 0.02 in agree- 
ment with the fluid turbulence case (Eq. (6) and Fig. 
5). Then, using Fig. 5 we can estimate the Taylor- 
microscale Reynolds number of the primordial turbulence 
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FIG. 6; As in Figs. 1,2 but for the CMB fluctuations in the 
3 year WMAP (Fig. 6a) and for the corresponding Gaussian 
spectral simulation map (Fig. 6b) . 
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FIG. 7: As in Fig. 3 but for the CMB fluctuations in the 3 
year WMAP (Fig. 7a) and for the corresponding Gaussian 
spectral simulation map (Fig. 7b). 



in the last scattering surface (with logarithmic accuracy) 
Rex - 10^. 

It should be noted, that critical value (for transi- 
tion from laminar to turbulent motion) is Rex ~ 40 
1^. Therefore, one can conclude that motion of the 
baryon-photon fluid at the recombination time is rather 
close to the critical state. This can be considered as an 
indication that the recombination process itself might 
be a source of the turbulence, which modulates the 
CMB. Indeed, the strong photon viscosity is known as 



the main reason for presumably laminar motion of the 
baryon-photon fluid just before recombination time. 
Therefore, if the recombination starts from appearance 
of the small localized transparent spots, then activity 
related to these spots can result in turbulent motion 
of the entire bulk of the baryon-photon fluid at the 
recombination time. 

I thank K.R. Sreenivasan for inspiring cooperation. I 
also thank D. Donzis, C.H. Gibson and participants of the 
workshop "Nonlinear cosmology: turbulence and fields" 
for discussions. 
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FIG. 8: The vortex ring and its velocity. 
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Appendix A 

Thin vortex tubes (or filaments) are the most promi- 
nent hydrodynamical elements of turbulent flows 
The filaments are unstable in 3-dimensional space. In 
particular, a straight line-vortex can readily develop a 
kink propagating along the filament with a constant 
speed. To estimate the velocity of propagation of such a 
kink let us first recall the properties of a ring vortex [29| . 
Its speed v is related to its diameter A and strength F 
through 



— In ("A 



(Al) 



where f] is the radius of the core of the ring and A/2?7 ^ 1 
(see figure 8). 

If, for instance, a straight line-vortex develops a kink 
with a radius of curvature A/2, then self-induction gen- 
erates a velocity perpendicular to the plane of the kink. 
This velocity can be also calculated using (A.l). One can 
guess that in a turbulent environment, the most unstable 
mode of a vortex tube with a thin core of length L (in- 
tegral scale) and radius r/ (Kolmogorov or viscous scale), 
will be of the order A: Taylor- microscale Then, the 
characteristic scale of velocity of the mode with the space 
scale A can be estimated with help of equation (A.l). 
Noting that the Taylor-microscale Reynolds number is 
defined as 0| 



Rex 



{A.2) 



where vq is the root-mean-square value of a component 
of velocity. It is clear that the velocity that is more rel- 
evant (at least for the processes related to the vortex 
instabilities) for the space scale A is not vq but v given 
by (A.l). Therefore, corresponding effective Reynolds 
number should be obtained by the renormalization of the 
characteristic velocity in (A.2) |^, as 



vX 



— In — 



^ ly 2ttv \2t] 
It can be readily shown from the definition that 



A 



(A.3) 



— — const Re 
where const = 15^/^ ~ 2. Hence 



1/2 



Re 



eff 



The strength F can be estimated as 
F ~ 27rt;,j77 



(^.4) 



(^.5) 



(^.6) 



where Vr^ = v/rj is the velocity scale for the Kolmogorov 
(or viscous) space scale 77 [igj. Substituting (A. 6) into 



(A. 5) we obtain 



i?e^^^ ~ ln(i?eA). 



Thus, for turbulence processes determined by the vortex 
instabilities the relevant dimensionless characteristic is 
ln(i?eA) rather than Re\ (cf Eq. (7) and Fig. 5). 

Appendix B 

In this Appendix we will show that modulation of the 
CMB by the primordial turbulence at the last scatter- 
ing surface can exhibit Gaussian properties much more 
strong than those of the 'normal' fluid turbulence, while 
preserving the turbulence properties intact. 

We always observe the CMB photons emitted from the 
'boundary' (the last scattering surface). Therefore, defi- 
nition of the 'boundary conditions' for the dynamic uni- 
verse is crucial in the case of the CMB radiation. Nat- 
urally this is a very non-trivial task to define boundary 
conditions for the universe (horizon). In the numerical 
simulations of the homogeneous turbulence dynamics the 
periodic boundary conditions are often used. The peri- 
odic boundary conditions seem to be also a natural sug- 
gestion for the universe. 

When one considers modulation of the CMB radiation 
by the Doppler effect only the relative difference in veloc- 
ity between the observer and the source, separated by a 
distance r along the line of sight: Au^ = u{x + r) — u{x), 
needs to be considered. Let us consider the moments of 
the velocity difference: 



Spir) = {[uix) - u{x + r)]P) 



(B.l) 
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FIG. 9: Moments Sp versus r / Li,ox obtained in a direct 
numerical simulation of fluid turbulence on a 2-n cube at 
f28^ resolution Re^ « 55 [s^l- Upper curves correspond to 
p = 2, 4, 6, 8, f 0, i2, the curves with symbols correspond to 
p = f , 3, 5, 7, 9, f f . For r/Lbox > 0.5 the picture is a mirror 
reflection of that shown in the figure. 



As an example, Fig. 9 shows results of a direct numeri- 
cal simulation [SOl of fluid turbulence with the periodical 
boundary conditions on a 2tt cube (size Lbox) at 128'^ 
resolution Re\ « 55. The curves in the upper part of the 
figure correspond to the even-order moments (B.l), while 
the curves with symbols correspond to the odd-order mo- 
ments. For r/I/boa: > 0.5 the picture is a mirror reflection 
of that shown in the Fig. 9. One can see that the Gaus- 
sian conditions: Sp ~ 0, for the odd-order moments are 
satisfied at the observer point r/Ljjox = 0.5. Actually, 
this is a direct consequence of the mirror symmetry pro- 
vided by the periodic boundary conditions. Using the 
dynamical equations of the baryon-photon fluid it can 
be shown that the even-order moments are also consis- 
tent with the Gaussian conditions at the observer point 
r/Li,ox = 0.5. 



